Interaction-induced excited-band condensate in a double-well optical lattice 
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We show theoretically that interaction effects in a double- well optical lattice can induce conden- 
sates in an excited band. For a symmetric double- well lattice, bosons condense into the bottom 
of the excited band at the edge of the Brillouin Zone if the chemical potential is above a critical 
value. For an asymmetric lattice, a condensate with zero momentum is automatically induced in the 
excited band by the condensate in the lowest band. This is due to a combined effect of interaction 
and lattice potential, which reduces the band gap and breaks the inversion symmetry. Our work 
can be generalized to a superlattice composed of multiple-well potentials at each lattice site, where 
condensates can be induced in even higher bands. 



Interaction between different orbitals lies at the heart 
of many challenging condensed matter phenomena such 
as heavy fermions and multi-band superconductors. 
Inspired by experimental developments in cold atom 
physics in the lowest band of optical lattices, physicists 
have begun to explore inter-band phenomena in optical 
lattices[T]-[6]. It is hoped that the precise controllability 
of optical lattices can provide unique opportunities to re- 
veal the interplay between multi-band physics and inter- 
action effects in strongly correlated systems. Although 
the specific results obtained in the multi-band optical 
lattice phenomena in these papers pertain to bosons, the 
generic considerations may also apply to fermions. 

An important question regarding the multi-band 
physics in optical lattices is how to populate the atoms 
in excited bands. Bosons have been excited to the p- 
band of a square lattice, where they establish coherence 
along specific directions [7 . Very recently, a novel tech- 
nique was used to populate bosons in the excited bands 
in a bipartite optical lattice [8]. p- and /- band con- 
densates were produced with life time up to a few tens 
of ms. These exciting experimental developments [71 [8] 
should enable the study of several interesting phenom- 
ena in higher bands of optical lattices recently predicted 
in the literature [9til3 . 

We emphasize that the condensates in excited bands 
obtained in current experiments are not the ground state 
of the system, and they eventually decay. It is then de- 
sirable to have a scheme for producing condensates in 
excited bands that have much longer life time. In this 
paper, we propose such a scheme using optical double- 
well lattices to produce excited-band condensates on top 
of an ordinary lowest-band condensate. The main results 
of this work are summarized as following. 

(I). In a symmetric double- well lattice, a condensate 
forms at the edge of the Brillouin Zone(BZ) in the lowest 
excited band with a finite momentum k, if the chemical 
potential is above a critical value that is achievable when 
the intra-band interaction is larger than the inter-band 
interaction. 



(II) . In an asymmetric double- well lattice where the 
inversion symmetery is broken, a stable condensate auto- 
matically forms at the top of the excited band with k = 
in the presence of a lowest band condensate. Unlike the 
case in (I), there is no threshold for this condensate to 
emerge and it can be induced by any finite condensate in 
the lowest band. 

(III) . Both the condensates discussed in (I) and (II) 
are the ground state properties of the system. They are 
stable and do not suffer a short life time problem. 

(IV) . Our studies on the double- well lattice can be gen- 
eralized to other multiple-well (more than two) super- 
lattices, for which even higher band condensates can be 
produced. 

(A) System and Hamiltonian: A double-well lat- 
tice can be produced by adding a short lattice of twice 
the periodicity along the x direction to the standard cubic 
lattice [14]. Physics discussed here can be easily general- 
ized to other lattice geometries [1 5 . The general form of 
the potential for the case we consider here can be written 
as 

where i = x^y, d is the lattice spacing, is the ampli- 
tude of the long lattice along different directions and Vs 
is the the amplitude of the short lattice along the x di- 
rection, cp characterizes the relative position between the 
two lattices. The symmetric case corresponds to cp = 
or 27r. 

When Vs is large enough, the lowest two bands are 
well separated from all the other ones. The many-body 
Hamiltonian in momentum space can be written as 

crk cr^ki 

(2) 

where a = g^e are the indices for the lowest two 
bands, e^rk = ^ak — ^crk is the single particle spectrum 
of each band and k is defined in the first BZ. U^^Z'^Z^V^ = 

kik2k3k4 
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FIG. 1: (Color online) (A) and (B) are schematic of a symmet- 
ric and an asymmetric double-well lattice respectively. The 
Bloch wave functions of the lowest two bands at k = are 
shown in red(bottom) and blue(top). (C) and (D) represent 
the pairing scattering = U^^_f^QQa,nd the interaction in- 
duced inter-band hybridization Gq* = 2^qqqq , which creates 
two particles at (k, — k) and one particle at k = in the ex- 
cited band respectively. The blue cloud represents the lowest- 
band condensate, and the holes represent the particles anni- 
hilated. 

^ / dR**^k, (R)* J.k. (R)*.3k3 (R)*.4k. (R)<5(ki + 

k2 — ks — k4 = K), K = 27r/(i(/i, ^2, ^3) is the reciprocal 
lattice vector, ^crk(R) = e*^ -^iicrk(R) is the Bloch wave 
function with crystal momentum k at band a, is the 
scattering length and M is the mass of the bosons. 

To be concrete, we focus on the Gross-Pitaevski 
regime, i.e., we assume there is a condensate at the bot- 
tom of the lowest band (a^o) = i^g- Because of the 
finite value of the gap, the higher band condensate is 
usually much weaker than the one in the lowest band. 
To the zeroth order, the magnitude of t/jg can be deter- 
mined by = (^^2^)^^^ where Ug = U^^^^ . For 
the excited band, there are many terms involved. Up to 
the order of ?/^^ and ?/^^, the Hamiltonian for the excited 

band can be written as He = X^^, Cgk^ek^ek -\- U -\- U' -\- 

Sk^o(^k(^ekS-k + Sk^e-k) + C.c) , 

U = ^IY^ (Wkaektte-k + C.c) , (3) 
k 

where eek = = U^^^l^, = U^^-i. 

and 

C/'=V^(GS4o + Goaeo), (4) 

where Gq = 2Z//qqqq . Because of momentum conservation, 
only contains a^Q or a^o. These leading order terms 
are shown in Fig. (IT]). The effect of higher order terms 
for determining the exact value of the excited-band con- 
densate amplitude will be discussed in section (B) and 
(C). 

(B) Symmetric lattices: Due to different parity 
of the wave function Ugo{R) and Ueo{R) in a symmet- 
ric double- well lattice, Gq vanishes, similar to the dis- 
cussion in reference [6 . We also find numerically that 



is generally one or two orders of magnitude smaller 
than for a broad parameter regions of lattice depth, 
i.e., Vl, Vs ~ 2 — IOEr. Therefore we ignore Wl, 
in the following discussions. Under this situation, 
is exactly the same as in the usual Bogoliubov the- 
ory for the excited band, i.e., = Xlk ^ek^ek^ek + 
Xlk (^k^ek^e-k + c.c). The grouud statc wave func- 
tion can be written as \G)p = Hk ^ «k '^ek^e-k |q^^ 



where = 




— 1, and Ee\c = 



^{U - M + 4V^Kokk)' - 4|I^kPV^^[16 . There is a crit- 
ical value /ic that can be written as 

^ic = + (4Z^ook?k. - 2|W^k, (5) 

where ki = 7r/(i(±l, 0, 0) is the location of the bottom of 
the excited band. If /i < /ic 7 the single-particle excitation 
in the excited band is gapped, i.e., Mm{£^ek} > as 
shown in Fig.([2[v). As a result, (bl^beji') ^ e"l^"^'l/^ 
decays exponentially, where 6eR = / dkctei^c^^^ is the 
annihilation operator in real space and the correlation 
length X ^ {1^ ~ Mc)~^^^- Though there is no single- 
particle condensate in the excited band, it can be shown 
that (^eR^eR) = — / (ik/^k/o^k ^ i^li i-e., there is an 
off-diagonal long range order in the reduced two-particle 
density matrix for the excited band. Therefore we refer 
to it as a pair-condensate in the excited band. 

On the other hand, by substituting /i = /i^ using 
Eq.([5| to the general expression for £^ek in the above 
paragraph, it is straightforward to show that E'gk = 
2|^ki|V^^\/^ek - ^eki if fi = fic is satisfied. When k 
ki, £^ek and the dispersion becomes linear at the 
edge of BZ as shown in Fig.(|2^). The vanishing gap indi- 
cates that a condensate emerges with a finite momentum 
ki. The condensate amplitude ip^^^ = (agki) can be cal- 
culated by minimizing the energy E = A\ipej^\'^ -\- B\ipej^\^ ^ 
where A = (U^ - M + ^i^Pooi'.u.) , B = and 
IV^cttI = —A/2B. While excited-band condensates have 
also been found theoretically in a standard optical lattice 
before O [4j, a double- well lattice has the advantage that 
a very large scattering length of atoms is not required 
since the band gap is significantly reduced. For example, 
we found that for Vl = ^Er, Vs = 2Er, Vz = I^Er and 
as/d = 0.04, V^^o/V^^ is about 3%. 

As /i is determined by Ug and t/jg through the relation 
M = ^^^0 + 2/7^?/^^, it is worthwhile to discuss under what 
conditions /i > /ic can be satisfied. For this purpose, we 
define Eg = - ^go > 0. Note W^,, = U^^'i.,^ ^ek, 
can be rewritten as 

Ecu, = + 2Eg{4c - 2)Ug^l^l + {I2c^ - 16c + 4)l7g>4, 

(6) 

where c = U^^^^^^/Ug. A straightforward analysis shows 
that if c < 1, £^ki vanishes at a sufficiently large value ipg. 
If c > 1, the single-particle excitation gap for the excited 
band will not vanish, no matter how large tijg is. The 
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FIG. 2: (Color online) The excitation spectrum for the 
ground (bottom from left) and the first excited band(top from 
left) of a symmetric double- well lattice for (A) /j < /Jc and (B) 



fic, where Vl 



VI = 15^H, VE = AEr Vs = IEr and 



as/d — 0.06. = 5 at /ic • Dashed lines in (A) represent the 
spectrum for non-interacting case. The ground band is also 
treated by Bogoliubov theory. 



reason is that increasing ipg has two competing effects 
on £^ek- It first increases the chemical potential (i ^ Ug 
which favors reducing the gap. On the other hand, it 
also increase the intra-band density-density interaction 
^ook-k- -"-^ other words, it enhances the single-particle 
gap and increase the threshold for the chemical potential 
to close the gap. Which effect wins depends on the ra- 
tio between Ug and Uqq^- Though we focus on the case 
where there is a large condensate in the lowest band in 
this work, our discussion can be easily generalized to the 
case where atoms in the lowest band form a Mott insula- 
tor instead of a condensate (see Supplementary Materials 
at [] for the discussions on Mott insulators). 

(C) Asymmetric lattices: Because the inversion 
symmetry is broken in asymmetric lattices, Gq becomes 
finite. It thus leads to an inter-band scattering process 
~ a\(^a\(^aaoCLeQ + c.c. We refer to this as an interaction 



induced inter-hand hybridization. In the mean- field de- 
scription, it gives rise to the terms in Eq.Q. 

The significance of U' is that it leads to a finite con- 
dense V^eo = (^eo) provided that ipg ^{). Unlike symmet- 
ric lattices where only the pairing term '0^(I^k^ik^i_k + 
c.c) is present, there is no threshold in the chemical po- 
tential for an asymmetric lattice to produce a condensate 
in the excited band. The only thing that matters is the 
ratio between ego — e^o and Go, which determines the 
amplitude of ?/^eO- An asymmetric double- well lattice has 
two unique advantages. First, it reduces significantly the 
energy difference between the lowest two bands Cgo — e^o , 
similar to what was found in [6 . Second, it provides 
flexibility in tuning the value of Gq. As a result, it is 
possible to achieve a reasonably large condensate in the 
lowest excited band, as we show below. 

To study how the condensate forms at k = in the 
excited band, we separate the following piece from the 
Hamiltonian in Eq.([3|, 

i^eO = /CeO + VeO + (V^o^OaL^L + V^aGoaJo + ^.c) (7) 
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FIG. 3: (Color online) (A) For an asymmetric double- well 
lattice, '0eo and as a function of '0pGo/(eeo — e^o), where 
V^oV^p/(eeO - Ig^) = 0.2 and [/eo/(eeO - e^o) = 0.4. (B) 
[/^(red solid), [/eo(blue dashed), 4^qqqq (purple dash-dotted) 
and Wo (Brown dotted) as a function of Lp that is defined 
in Eq.([T]) to characterize the relative position of the long 
and short lattice potential, where V^/Er — VI /Er — 15, 
V£/Er = 2, Vs/Er = 6, and as/d = 0.06. (C) Go (red solid) 
and ^eo — ^^o(blue dashed) as a function of Lp. The param- 
eters are the same as (B). (D) V^eo/V^^ for the same lattice 
configuration as (B) and (C). 



where /Cgo = eeoa^Qaeo and Veo = ZYoooo^L^L^eOtteO. 
Without losing generality, a repulsion Vgo has been added 
to the Hamiltonian. We have also set both iXgO and Ug^ to 
be real. There are two different driving terms for the ex- 
cited band condensate. The quadratic term ttgo^eo f^^vors 
the formation of a pair-condensate. The linear term a^Q, 
however, induces a single-particle condensate. The com- 
petition between them leads to an interesting evolution 
from a pair-condensate to a single-particle condensate in 
the excited band with increasing Gq. 

The ground state for the Hamiltonian in Eq.([7|) is writ- 
ten as \G)s = Zln^^"^!^)' where the coefficient 
can be calculated exactly. The order parameters for 
the single-particle and pair condensate are defined by 
V^eO = (^eo) and ip'^Q = (ago^eo) respectively. An illumi- 
nating example to demonstrate how i/jeo and V^^q change 
as a function of tpgGo is shown in Fig.(|3jA). As seen from 
Eq.Q, the amplitude of i/jeo depends on the product of 
t/jg and Gq. On the other hand, it also depends on the 
value of Ceo — e^o- The larger the latter is, the smaller 
the amplitude of ipeo is. Thus, we have chosen the di- 
mensionless number ^pgGo/{ceo — ^go) as the x-axis in 
Fig.([3|^). When Go = 0, V^eO = 0, we have a pure pair- 
condensate. The negative value of i/j'^q simply means that 
this pair condensate has a relative phase tt with respect 
to t/jg^ since Wo > in our case. With increasing Go, 
IV^eol increase while IV^^qI decreases. In the large GoV^^ 
limit, ipeo I V^eO 1, where ngo = (^L^eo)- This indi- 
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cates that the ground state for the Hamiltonian in Eq. (|3| 
becomes a coherent state. 

To concrete the above discussions for real cold atom ex- 
periments, we have carried out an exact band structure 
calculation for various double- well lattice configurations. 
The parameters in Eq.jTl for a typical lattice configura- 
tion are shown in Fig.(|3p,C). The condensate amplitude 
in the excited band is shown in Fig.([3)l)) for different val- 
ues of (p. We have set ip'^ = 2. It is clear that V^^q can 
be as large as a few to ten percent of ijj'^. Thus for total 
particle number TV ~ 10^ — 10^ in the weakly interact- 
ing region, the condensate particle number in the excited 
band can be of the order of 10^ — 10^. For the lattice 
configuration as shown in Fig(|3^-D), the maximum of 
T/^gQ occurs at = TT, where Go also reaches its maximum 
as shown in Fig.(|3p). We emphasize that this interac- 
tion induced condensate comes from the combination of 
interaction effect and the lattice potential that breaks 
the inversion symmetry. Without reducing the value of 
^eo— ^^^0, the condensate amplitude would be rather small 
though a finite coupling in principle can exist if the in- 
version parity is broken for any reasons. This interaction 
induced condensate in the excited band is a property of 
the ground state, since the interaction naturally mixed 
the condensates in the lowest two bands. Thus, it does 
not suffer from the short life time problem of the usual 
single excited band condensate. 

(D) Probing condensates in different bands: The 
crystal momentum distributions from different bands 
overlap with each other in the standard Time-Of- 
Flight(TOF) expansion. In a ballistic expansion, the in- 
terference pattern represents the original momentum dis- 
tribution (nk) = / ^^Ri^^R2e'^'^^'"^'^(^K^^R2), where 
nk is the density operator in the momentum space 
and Byi is the field operator in the real space. Ex- 
panding Byi in the basis of Bloch wave functions, i.e., 
Byi = X^crk ^crk(R)^crk, ouc immediately sees that 
(^k) = Eaa'pqCp(k)0a'q(k)(atpa^,q), whcrc 0aq(k) = 

J (iRe~*^ "^^crq(R')^^k-q=K IS the Fourier transform of 
the Bloch wave function. Therefore, the population at 
the reciprocal lattice vectors k/ = /(27r/(i, 0, 0), where / 
is an integer, can be written as 

= Yl Co(k0^a'o(kO(4o^a'o). (8) 

a=g,e 

Thus one can first numerically obtain the Bloch wave 
functions ^crk(R) and its Fourier transforma ^crq(k) for 
the double- well lattice, and then use them to fit the TOF 
images. From the amplitude of a number of peaks in the 
interference pattern, one is able to extract (a^Qacr'o), i-e., 
the condensate amplitude at k = of different bands. 
For the condensate at the edge of BZ, additional peaks 
located at k = 7r/dx-\-l'K emerge as the unique feature of 
this condensate with a finite crystal momentum. Another 
method is to apply the band mapping technique [17|, 



which can directly map out the crystal momentum dis- 
tribution where condensates of different bands can be 
directly measured. 

(E) Condensates in a superlattice: Finally, we 
discuss the general principle of producing condensates in 
excited bands of a superlattice composed by multiple- well 
potentials in each single site. (1) A simple way to create 
a superlattice containing n wells at each site is to use 
two lasers with wave lengths A and nX. The lowest few 
bands may be well separated from other higher bands. 
(2) When this multiple-well potential respects inversion 
symmetry, the lowest band condensate automatically in- 
duces a condensate in an excited band of even parity at 
k = 0, because of the presence of the interaction term 
U{ = G^a^gQCL^gQagoae'o + c.c in the Hamiltonian, where 
C ^ 0. For odd-parity bands, and U{ vanish. How- 
ever, U2 = ^'^^Qa^Qae'-ki^e'ki + c.c is always present, 
and if the chemical potential is large enough to overcome 
the single-particle excitation gap for those bands, con- 
densates form at the edge of BZ with a finite momentum 
ki . When the inversion symmetry is broken for the asym- 
metric multiple- well potentials, G' becomes finite for all 
the excited bands, and condensates in the lowest few ex- 
cited bands can be induced. 

In summary, we have pointed out that in a double- 
well lattice, a stable condensate in the lowest excited 
band can be induced by inter-band interactions. It ei- 
ther comes from large enough interaction that overcomes 
the single-particle excitation gap for the excited band, or 
a combination of inter-band interaction effect and the lat- 
tice potential that breaks the inversion symmetry. Our 
proposal can be generalized for producing condensates 
in the lowest few excited bands of a superlattice. Re- 
cently, Sengstock's group has reported an observation 
of interaction-induced condensate in an excited band of 
a triangular double- well lattice [18 . This discovery di- 
rectly proves the validity and feasibility of our proposals. 
This work is supported by JQI-NSF-PFC, ARO-DARPA- 
OLE, ARO-MURI and AFOSR-MURI. 
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